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ABSTRACT 
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Aims. We have performed optical polarimetric observations of the SN 2006aj associated to the X-ray flash (XRF) of February 18, 2006, 
\Q XRF 060218 that provide information on its expansion geometry. 

f"*) . Methods. The data were acquired in the iJ-band with the 0.7m telescope of Crimea, 2.5m Nordic Optical Telescope and the 2.2m of Calar Alto. 
Results. We report the detection of linear polarization between 3 and 39 days after the gamma-ray event (/ - to). This represents the first 
polarization detection of a Ic supernova (SN) associated to an XRF. Our data exhibit a degree of linear polarization (P) around P ~ 4% at 
t - to ~ 3 - 5 days, followed by a constant polarization phase with P ~ 1.4% at 13.7 < / - to < 39 days. Our data suggest a decay in P, and 
more interestingly, show a position angle {6) rotation of ~ 100° comparing data taken before and after the R-hanA lightcurve peak. 
Conclusions. The reported polarization measurements can be explained by the evolution of an asymmetric SN expansion. We discuss on several 
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^ ■ ingredients that could account for the observed 6 rotation. 

C3 • 

^ ■ Key words. Gamma rays: bursts - Supernovae: general - Techniques: polarimetric 

d ' 1. Introduction some counterexamples exist (Leonard et al. 2006). The very 

few polarimetric studies of Ic SNe related to GRBs do not al- 
The Swift satellite detected a long-soft gamma-ray burst (GRB) 1qw ye{ tQ infer ^ global polarimetric differences with nor- 
at k) = Feb. 18.149 2006 UT, (Cusumano et al. 2006), showing ma] fc SNe However5 as we report in this Let ter, some normal 
both a low redshift (z = 0.033, Mirabal & Halpern 2006) and fc SNg ^ SN 2002ap (Mazzali et a] 2QQ2) ^ SQme sim _ 
a peak energy (Campana et al. 2006), typical of X-ray flashes ^ SN 2006aj . polarization leyels above „ 1% andi 

(XRF 060218 hereafter). Further observations associated it to a • . , , . ... ... , • .. , 

v ' . . _ more interestingly, long-term rotations in the polarization angle 



Ic-type supernova, SN 2006aj (Masetti et al. 2006 1 Soderberg et (Kawabata et al 2002 1 
al.|2006| Plan et al.|2006J. This finding strengthened even more Most of the S N intrinsic polarization detections are at- 
the already solid lon g-durat ion GRB/Supernova link (Hjorth et tribufgd to the break of the symmetry around line of sight, due 
alHOOll Stanek et al.| 2003 |l. to a Thompson photon scattering through an aspherical SN ex- 
Since 1999 (Covino et al.|1999) several GRB afterglows pan sion (Kasen et al. 2003 1. This polarization is expected to 
have shown op tical p olarization at a level of -1-3% (see decay as ~ r 2 , as the atmosphere expands and t he opti cal depth 
Gorosabel et al. [2004J and references therein). Optical super- (and hence the polarization) falls (Leonard et al. l2006ll . In most 
novae (SNe), including typical Ic-types not related to GRBs, cases fh e afterglow is so strong that overshines the SN polar- 
do not tend to show optical polarization above ~ 1%. However, ization component (for instance in GRB 030329; Greiner et al. 

2003). In this Letter we report the first polarization detection 



to date of a SN associated to an XRF, not coming from the 



Send offprint requests to: J. Gorosabel, e-mail: jgu@iaa.es 

* Based on observations made with the Nordic Optical Telescope, 

operated on the island of La Palma, in the Spanish Observatorio del GRB ^ecte forward shock emission (Fan et al. 2006, Soderberg 

Roque de los Muchachos of the Instituto de Astrofisica de Canarias, et al - ESS- Sect. EJ reports the observations earned out for 

and on observations done at the Centro Astronomico Hispano Aleman SN 2006aj . Sect. presents our results and Sect. H discusses 

(CAHA) at Calar alto, operated jointly by the Max-Planck Institut fur the implications of our measurements. Sect. summarises the 

Astronomie and the Instituto de Astrofisica de Andalucfa (CSIC). final conclusions. 
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2. Observations and data analysis 

Table [^displays the log of our observations. All the measure- 
ments reported in this paper were earned out in the R -band. The 
observations were performed, chronologically ordered, with 
the 0.7m telescope of Crimea (AZT-8), the 2.5m Nordic Optical 
Telescope (NOT) and the 2.2m of Calar Alto (CAHA). 

The AZT-8 observations were done with two Savart plates; 
swapping them the observer can obtain either the Stokes Q (two 
images are split diagonally) or the U parameter (split horizon- 
tally). The field of view (FoV) covered by the AZT-8 (+ST-7 
CCD) is 8.1' x 5.4' and the pixel scale is 1.3"/pix. Since the 
GRB field is substantially crowded in the AZT-8 images, we 
performed PSF photometry (Stetson fl987> . 

The NOT observations were based on ALFOSC equipped 
with the FAPOL unit. In order to obtain Q and U with 
ALFOSC, the GRB field was imaged through a calcite and a 
1/2 wave plate. Four images of the GRB field were acquired, 
rotating the 1/2 wave plate at 0°, 22.5°, 45.0° and 67.5°. The 
calcite plates of FAPOL produce a vignetted field of about 140" 
in diameter with a pixel scale of 0.19"/pix. The CAHA obser- 
vations were carried with CAFOS. The CAFOS polarization 
unit is similar to FAPOL (is based on a Wollaston prism in- 
stead of a calcite plate), but uses a strip mask on the focal plane 
to avoid accidental overlapping on the CCD. The total FoV in 
CAFOS is composed by 14 strips of 9' x 18" each, contain- 
ing stars enough for a satisfactory interstellar medium (ISM) 
polarization correction. For the NOT and CAHA data the de- 
termination of P and 9 was carried out fitting the S (6) function 
with the corresponding ISM normalization factor (di Serego 
Alighieri 1997, see Fig. [3- For both data sets aperture pho- 
tometry was performed, ranging the aperture radii from 1 pixel 
to 2.5 times the FWHM (Full Width at Half Maximum). After 
checking that the magnitudes were independent of the apertures 
used, the ones yielding minimum errors were adopted. 

The Galactic ISM polarization correction was done using 
9 (13) field stars on the AZT-8 (CAHA) images. Due to the 
reduced ALFOSC FoV only one unsaturated bright field star 
was available in the NOT data. Thus for the NOT images 
zero (G191B2B) and high (HD25443) polarization standards 
allowed us to infer the SN 2006aj Stokes parameters, and for 
completeness, also the ones of the unique unsaturated field star 
(which agreed with the ISM law by Serkowski et al. U975> . 

3. Results 

Table[2shows the position angle (0) and the degree of linear po- 
larization (P) of the afterglow once they were corrected by the 
ISM polarization and also by the fact that P is a positive quan- 
tity (multiplying P by y/l - (crp/P) 2 ; see di Serego Alighieri 
I1997> . The AZT-8 results displayed in Table [0 supersede the 
ones reported previously (Larionov & Larionova 2006 1. The 
four AZT-8 data sets yield P marginal detections, reaching a 
3cr confidence level at t - to ~ 3.6 days. The apparently large 
P variations (from P=2.43 to 4.53) in the first four nights are 
within P errors, so we they are likely statistical fluctuations. 
However, we note that for the four epochs SN 2006aj is sys- 
tematically placed at the same Stokes plane area, once the ISM 
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Fig. 1. Upper panel: The Q, U plane of the stars in the field 
(circles) of XRF 060218 (square) imaged from CAHA on 
March 8.845 UT (t - to = 18.7 days). The size of each cir- 
cle is proportional to the distance of the corresponding star to 
the center of the CCD. This is used to check that the instru- 
mental polarization, if present, does not vary severely across 
the FoV. Lower panel: The satisfactory fit of the S ((f)) function 
(x 2 /d.o.f = 1.38). As shown the XRF 060218/SN 2006aj is 
clearly polarized. An aperture radius = 1 FWHM was assumed. 



correction is included. In order to reinforce the P detection, the 
four AZT-8 observing epochs were co-added and PSF photom- 
etry performed. The fifth line of TableQ] shows the result when 
the four AZT-8 epochs are combined. The estimate of the P 
and 6 errors in the AZT-8 images were obtained modelling with 
a Montecarlo method the flux error distribution of SN 2006aj 
and the 9 field stars used for the ISM correction. Fig. [2] shows 
the Stokes plane when all the AZT-8 data are combined. The 
cloud of dots representing the SN 2006aj (Q, U) distribution is 
clearly off the center defined by the field stars. The P value in- 
ferred from the combined image is 3.40 + 1 .09%, detected at a 
3. lcr level (including the yjl — (crp/P) 2 term). 

Both the NOT and first CAHA data show conclusive P de- 
tections with confidence levels of 10.3<x and 6.7<x respectively. 
The 6 values of both epochs agree with a stable angle around 
50°, but interestingly it appears rotated by ~ 100° with respect 
to the AZT-8 images. Our data points reject a stationary P value 
at a 3<x level, so we conclude that P likely fades at 5.6 < t— to < 
13.7 days. Such a P decay resembles the late evolution of the 
asymmetric SN 2004dj (Leonard et al. 2006 ), which showed an 
abrupt P peak at t - to ~ 90 days due to the sharp appearance 



of the inner ejecta, followed by a late 



decay. 



We point out that the result of our /?-band polarimetry is 
a weighted mean of the continuum and the spectral features 
(in principle each one with a different polarization) falling in- 
side the wavelength range of the filter transmittance (centred 
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Table 1. Log of the observations. For the NOT and CAHA the observation are cycles of 4 images, (4 x T exp ) with the plate 
rotator at 0, 22.5, 45, 67.5 degrees. The R-bwA magnitudes are calibrated using the standard stars given by Hicken et al. (2006 1. 



Days from GRB 


Telescope 


Exposure Time 


Seeing 


Magnitude 


P 


e 


t- 


h 


(+Instrument) 


seconds 




tf-band 


% 


degrees 


2.569954 - 


2.597558 


0.7m AZT-8(+ST-7) 


2 x (10 x 90) 


3.0 


17.99 + 0.03 


2.61 ± 1.80 


157.82 ± 16.47 


3.579178 - 


3.600764 


0.7m AZT-8(+ST-7) 


2 x (10x90) 


3.8 


17.88 + 0.02 


4.53 ± 1.52 


143.64 + 7.18 


4.555046 - 


4.599144 


0.7m AZT-8(+ST-7) 


2 x (20 x 90) 


3.3 


17.68 + 0.02 


4.06 ± 1.63 


143.86 + 7.98 


5.565510- 


5.612326 


0.7m AZT-8(+ST-7) 


2 x (20 x 90) 


4.2 


17.52 + 0.02 


2.43 ± 0.89 


150.09 ± 10.02 


2.569954 - 


5.612326 


0.7m AZT-8(+ST-7) 


2 x (60 x 90) 


3.5 


17.71 ± 0.02 


3.40 + 1.09 


147.56 ± 7.06* 


13.721516 - 


13.755185 


2.5m NOT(+ALFOSC) 


4x700 


1.0 


17.19 + 0.04 


1.44 ±0.14 


49.58 ± 2.73 


18.673955 - 


18.719135 


2.2m CAHA(+CAFOS) 


4x900 


1.6 


17.38 + 0.03 


1.40 + 0.21 


51.62 ± 4.11 


38.679248 - 


38.707514 


2.2m CAHA(+CAFOS) 


4x400 


2.0 


18.44 + 0.03 


< 3.9 1 ' 





★ Determined combining the images of the above four observing epochs. 
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Fig. 2. The Stokes plane of 9 field stars and SN 2006aj on 
the combined AZT-8 image when the Galactic ISM correc- 
tion is applied. The clouds of dots were created altering with 
a Montecarlo method the fluxes of the objects according to 
their photometric errors. The cloud representing SN 2006aj is 
shifted with respect to the origin defined by the field stars, as- 
sumed to be intrinsically unpolarized in average. 



at 6500A). With no spectro-polarimetric data it is not possi- 
ble a priory to quantify if the measured polarization evolution 
is dominated by spectrally localized large polarization fluctu- 
ations. We note however that the SN 2006aj spectra around 
6500 A are smooth. Thus the 6 rotation measured at 5.5 < 
t - to < 13.7 days, is not easily explainable by colour evolu- 
tion, which occurs at wavelengths bluer than 5000A in that 
epoch range (Mazzali et al. 2006 1. 

4. Discussion 

The most natural framework to explain our polarization mea- 
surements is a non spherical SN expansion, consistent with the 
scenario proposed for XRF 0602 1 8 by recent studies (Pian et al. 
2006 1. It is interesting to note that other Ic SNe exhibiting non- 
spherical expansions have been linked to long-duration GRBs 
(Hoflich et al. [T9991 Further more there is growing evidence 
that an important fraction of Ic SNe are triggered by bipolar 
ejections (Granot & Ramirez-Ruiz 200?). It is widely accepted 
that non-zero SN polarization measurements demand some de- 



gree of expansion asphericity. If the projection of the SN along 
the line of sight is elliptical, a non canceled linear polarization 
should be present parallel to the minor axis (Kasen et al. 2003 1. 
An ellipsoidal expansion is fairly consistent with the measured 
P (see r 2 curve of Fig.|3}. However, a single ellipsoidal model 
can not account for the measured 8 rotation. 

The decay in P seems to be accompanied with a decay of 
the photospheric expansion velocity (notice the similarity be- 
tween our Fig. [3] lower panel and Fig. 2 of Pian et al. 2006 1. 
This possible correlation between the photospheric expansion 
speed and the polarization, might imply the existence of a fast 
(v > 20000 km s _1 , Pian et al. 2006 1 highly asymmetric ejecta 
responsible of the high polarization measured 2.6 - 5.6 days 
after the GRB. This rapid asymmetric ejecta might be aligned 
along the GRB jet. Several days later a slower (v < 15000 
km s -1 , Pian et al. 2006 1 SN bulk ejecta would circularice the 
geometry, which would explain the lower P of our CAHA and 
NOT measurements. We could speculate that this slower ejecta 
might be in a toroidal-shape equatorial expansion, so the geo- 
metric transition from the rapid to slow ejecta might naturally 
explain the detected rotation of ~ 100°. These results partially 
resemble the jet-like geometry proposed for SN 2002ap, which 
exhibited lower P values than SN 2006aj, but with a Q rotation 
similar to it (~ 70°; Kawabata et al. 2002 1. However, other in- 
gredients can also induce a rotation in 0: i) the host galaxy ISM, 
it) a non radial distribution of 56 Ni and Hi) a dusty circumstel- 
lar medium (CSM). These effects can be combined producing 
a complex picture. 

As discussed by Klose et al. ( 2004 1 the polarization mea- 
sured in afterglows are sensitive to the host galaxy ISM, so 
it might well affect our observations. Thus, the 9 stability at 
t - to > 13.7 days could be the result of the P induced by 
the host ISM, so the SN 2006aj intrinsic polarization (Psn) 
would dominate just before our NOT data point. According to 
the Serkowski's law E(B — V)> 0.15 would be required in the 
host line of sight to account for the P ~ 1 .4% NOT and CAHA 
measurements. This E(B - V) value agrees fairly well with the 
host extinction estimated by Modjaz et al. (2006, E(B - V) ~ 
0.05 -0.11), So llerma n et al. d20Tj6l E(B - V) ~ 0. 1 - 0.3), and 
Campana et al. J2006I E(B - V) = 0.2). We note however the 
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lower estimate by Guenther et al. (EUUfjl E(B - V) = 0.042 1 ), 
so the extinction although present, might not account totally for 
the late polarization. In this Framework the 9 change of ~ 100° 
would just obey to the chance alignment between the asymmet- 
ric expansion and the linear polarization induced by the host 
ISM. Assuming that our last P detection (? - to ~ 18.7 days) is 
purely due the host ISM, the SN 2006aj intrinsic polarization 
at t -t ~ 4.5 days would be higher (P SN = 4.78 + 1.17%, see 
Fig. flower panel). 

An asymmetric distribution of 56 Ni could alter the SN ra- 
diation field, specially after the lightcurve peak when the 56 Ni 
radioactive decay powers the optical lightcurve. The collapsar 
model predicts a subrelativistic disk wind of neutron and pro- 
tons, which after cooling would form 56 Ni (Woosley & Bloom 
20061. A picture to explain the 9 rotation of ~ 100° 



would 

require an equatorial deposition of synthesized 56 Ni, which 
would alter the 9 measured after the lightcurve peak. 

Another possible mechanism to explain the observed 9 
change could be delayed CSM scattered light. This process 
was proposed to interpret the SN 1987 A polarization lightcurve 



(Wang & Wheeler 119961 1. In this scenario the polarization 
would consist of two components, the one due to photons 
propagating directly from the (possibly aspherical) SN, and 
the second one from scattered light by a dusty CSM region. 
The epochs of our observations would push the CSM scatter- 
ing region further than ~ 2.3 light-days from the progenitor. 
However, this process should have originated a clear distortion 
of the spectra at A > 6000 A after the lightcurve peak, which 
is not observed (Modjaz et al. 2006 Pian et al. 2006 1. 
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Fig. 3. Our polarization results, once the the Galactic ISM cor- 
rection was considered. Upper panel: As shown 9 exhibited a 
rotation of ~ 100° at 5.6 < t - to < 13.7 days. Middle panel: 
P shows a likely decay reaching a constant level at a ~ 1.4% 
level. Lower panel: The intrinsic P$n evolution assuming that 
our last P detection is entirely originated by the host ISM. The 
dashed curves show the +lcr uncertainty region of the f 2 de- 
cay, as seen in SN 2004dj. 
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5. Conclusions 

We report the detection of 7?-band linear polarization for 
XRF 060218/SN 2006aj in observations done since t - t ~ 3 
to 39 days. These observations represent the first detection of 
optical polarization in an XRF. The measured high degree of 
polarization gives further credence to the XRF/GRB off-axis 
scenario, since a face-on SN would not likely produce any in- 
trinsic polarization. As Patat et al ( 2001 1 noted there is a degen- 
eracy between the viewing angle and the asymmetry level, so 
a detailed geometric picture based on only our data is not pos- 
sible. Our data show a P ~ 4% value in the first ~ 3 - 5 days 
after the GRB, followed by a constant P ~ 1.4% polarization 
phase at t - to > 13.7 days, suggesting a P decay. 9 shows a 
~ 100° rotation at 5.6 < t- to < 13.7 days. We propose that the 
P evolution can be explained by an asymmetric SN. We discuss 
several ingredients which separately, or combined, can explain 
the observed 9 rotation (in decreasing order of relevance); i) 
a highly asymmetric high-velocity SN ejection, followed by a 
less asymmetric SN bulk, ii) polarization induced by the host 
ISM, iii) asymmetric distribution of 56 Ni and iv) scattering by 
a dusty region placed a few light-days from the progenitor. 



1 See Sollerman et al. 1 2006 1 for a discussion on possible uncertain- 
ties and biases affecting this measurement. 
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